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ABSTRACT 

By integrating the relativistic hydrodynamic equations introduced by Taub we have determined the exact EQ- 
uiTemporal Surfaces (EQTSs) for the Gamma-Ray Burst (GRB) afterglows. These surfaces are c ompared an d 
contrasted to the on es obtained, using approximate methods, bv lPanaitescu & Meszarosl d!998h : ISaril (TT998): 
lC,ranotetalJ ( rf999i) . 

Subject headings: black hole physics — gamma rays: bursts — gamma rays: observations — gamma rays: 
theory — relativity 



1. INTRODUCTION 

The recent explanation of the observed luminosity in X- and 
7-ray energy bands in Gamma-Ray Bursts (GRBs), as well 
as the comprehension of their spectral properties, depends in 
an essential way on the determination of the EQuiTemporal 
Surfac es (EQTSs) in the afterglow era (iRuffini et al J 120021 
2003a). Here we compare and contrasts the exact determina- 
tio n of the EQTSs with the approx imate e xpress i ons presented 
bv iPanaitescu & Meszarosl il998l> : ISaril (119981) : iGranot et al.l 
(1999). 

A great deal of consensus has been reached concerning 
three basic issues in the description of GRB afterglows: 

a) Their origin is generally traced back to the interaction of 
an ultrarelativistic baryonic matter pulse with the InterStellar 
Medium (ISM). It is also agreed tha t the r elativistic hydrody- 
namic equations introduced by Abe iTaubl (TT948 ) are the cor- 
rect theoretical framework to describe such an interaction. 

b) The general definition of the EQTSs is also agreed upon by 
everyone. 

c) The necessity of determining the boundary conditions of 
the baryonic matter pulse in the early phases of the afterglow 
era is also generally recognized. 

We illustrate in the following sections the equations needed 
for the description of these three basic issues and also point 
out a major difference between our approach and the ones 
in the current literature concerning the solutions of the Taub 
equations. We identify in this difference the reason for the re- 
sults on the EQTSs presented in Fig.[3]and Fig.|4] Our treat- 
ment assumes the exact integration of the Taub equations. 

We recall that all the GRB observable quantities depend es- 
sentially on the EQTSs. Due care is therefore needed in their 
correct determination. 

2. THE THEORETICAL BACKGROUND 

2.1. The Taub relativistic hydrodynamic equations 

The relat ivistic hydro dynamic equations have been ex- 
pressed by ITaubl {l948) in the analysis of the relativis- 
tic Ra nkine-Hugoniot equations, see also lLandau &~ Lifsits 
( 119951) for the explicit comparison between the relativistic 
and nonrelativistic regimes. Although consensus has been 
reached on the equations to be used in the description of 
the GRB phenomenon, differences still exists on the va- 
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lidity of the approximations adopted in the solutions and 
their physical interpretation. The Taub equations have been 
used to describe the adiabatic optically thick expansion phase 
of the pulse generating the GRB (see e.g. MeszarosetalJ 



1993 



Bisnovatyi-Kogan & Murzina 1995; Ruffini et al. 199S 



2000). In the present context of the GRB afterglow con- 



sensus exists that these equations acquire the fo rm (see e. g. 



Blandford & McKeel1976tlpir^rIT999: Rufnnjj£id|2003b): 
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where E mt , 7 and M are respectively the internal energy, the 
Lorentz factor and the mass-energy of the expanding pulse, 
Tiism is the ISM number density which is assumed to be con- 
stant, m p is the proton mass, e is the emitted fraction of the 
energy developed in the collision with the ISM and M lsm 
is the amount of ISM mass swept up within the radius r. 
Msm = (4/3)7r(r 3 - r 3 )m ; ,ni sm , where r is the starting ra- 
dius of the shock front. In general, an additional condition 
is needed in order to determine e as a function of the radial 
coordinate. In the following, e is assumed to be constant and 
such an approximation appears to be correct in the GRB con- 
text. 

iBlandford & McKed iwffl) obtained a self-similar solution 
of the Taub equations in a so-called "ultrarelativistic" approx- 
imation with: 

7o > 7 > 1 , (2) 

where 7 is the initial value of the Lorentz gamma factor of 
the shock front. In the current literature, such an approxi- 
mation has been used in order to obtain a simple power-law 
relation between the Lorentz gamma factor and the radial co- 
ordinate of the expanding shock: 



7 oc r 



(3) 



where a = 3 corresponds to the fully radiative case (i.e. e = 1 
in Eqs.Q) and a = 3/2 to the fully adiabatic case (i.e. e = 
in Eqs.Q). 

In our approach we have exactly integrated the Taub equa- 
tions. For all the sources which we have analyzed, including 
the case of GRB 991216 which will be used as a prototype in 
this paper, we have always found j a ~ 200 - 300. Under these 
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conditions, we have found that no power-law solutions like 
Eq.Q hold in any finite region of the afterglow: at most we 
can define at every point an "effective" power-law index a e s 
which in the early phases of the afterglow is zero, then reaches 
a maximum value and then decreases all the way back down 
to zero in the latest phases of the afterglow. Such a maximum 
value of a e g is always found to satisfy a e ff < 3 in the fully 
radiative case and similarly a e ff < 3/2 in the fully adiabatic 
case (see Fig.^below). 

2.2. The definition of the EQTS 

The EQTSs are surfaces of revolution about the line of 
sight. The general expression for their profile, in the form 
$ = i?(r), corresponding to an arrival time t a o f the ph otons 
at the detector, can be obtained from (see e.g. R uffini et al.l 
I2003bl) : 

ct a = ct (r)-rcos$ + r*, (4) 

where r* is the initial size of the expanding source, ■& is the 
angle between the radial expansion velocity of a point on its 
surface and the line of sight, and t = t(r) is its equation of mo- 
tion, expressed in the laboratory frame, obtained by the inte- 
gration of Eqs.Q. From the definition of the Lorentz gamma 
factor 7 -2 = 1 - (dr/cdt) 2 , we have: 

ct = JZ[l-rHr)Y 1/2 dr, (5) 

where 7(r) comes from the integration of Eqs.Q. 

In the current literature, the initial size r* is usually ne- 
glected: this is certainly possible in the description of the 
latest phases of the afterglow, but not in the earliest ones 
dRuffini et alJl200 2). The big difference, however, is that in 
our approach we use the exact solution of Eqs.Q in the eval- 
uation of Eq.(|5j, while in the current literature the approxi- 
mate solution given in Eq.(|3} is genera lly used with a serie s 
of additional approximations (see also Ruf fini et al.l l2001a). 
All this leads to the results presented in Fig.[3]and Fig.|4] 

2.3. The boundary conditions 

Our model, like other ones in the current literature (see e.g. 
Mesz aros et al.119931) . is function of only two parameters: the 
total energy E to , of the electron-positron pairs originating the 
GRB phenomenon and their baryonic loading. In our model 
E tot coincides with the energy of th e dyadosphere E c { ya (see 
RuffinH U998l iPreoarata et aLlH998l) . The expansion of the 
photon and electron-positron pairs pulse has been computed 
in full details both by semi-analytic treatments and numer- 
ical sim ulations (the pair-electromagnetic [PEM] pulse, see 
Ruf fini et alJ ll999). In our model the baryonic load is due to 
the engulfment by the PEM pulse of the baryonic matter left 
over by the collapse of the progenitor star at a radius r ~ 10 10 
cm in the still optically th ick and adiabatic expansion of the 
pulse (Ruffini et al. 2000). The reaching of the transparency 
condition of this ele ctron-positron pair and baryonic matter 
is computed (see e.g. lBianco et alJ t2001). The determination 
of the only two free parameters is then made by fitting the 
average intensity of the rising part, of the peak emission and 
of the d ecreasing intensity part of the afterglow (Ruffi ni et al.l 
2001b). These parameters fixes uniquely the initial Lorentz 
gamma factor of the baryonic pulse, as well as its initial mass 
and the arrival time at the detector, t~ = (1 +z)t a , at which 
the afterglow begins. As usual, z is the cosmological red-shift 
of the source. In the specific case of GRB 991216 we have 
7o = 310.1, r Q = 1.943 x 10 14 cm, (f<') = 8.413 x 10" 2 s. 




Radial coordinate (r) (cm) 

FIG. 1. — The Lorentz gamma factors are plotted as a function of the 
radial coordinate both in the fully radiative and adiabatic cases for the exact 
(approximate) solution with solid (dashed) lines. 
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FIG. 2. — The radius r of the expanding pulse as a function of the arrival 
time computed at $ = is plotted both in the fully radiative and adiabatic 
cases for the exact (approximate) solution with solid (dashed) lines. 



3. THE TREATMENTS BY SARI AND PANAITESCU & MESZAROS 

In the current literature (see e.g. iPanaitescu & Meszaros! 
I1998I) the description of the earliest parts of the afterglow are 
neglected and special attention is given to later times when 
the effect of the deceleration becomes important and Eq.(|3} is 
assumed to apply. The following approximations have been 
adopted: 

a) Eq.Q is assumed to hold only during the so-called "de- 
celeration phase" when 7 < 7^ = (2/3)7 and r > r c {, where 
l = Jd(r/r d )~ a ; 

b) instead of Eq.0 the following approximate expression is 
used: 

ct = r d [l + (4a + 2)- 1 7 - 2 (r)]+^ [l + (l/2) 7 - 2 (r)] dr; (6) 

c) no treatment of 7 is given for r < r d . 

In order to avoid additional arbitrary factors, we compare 
the results by adopting in the above approximate equations 
the values of 70 and obtained from our exact solution. In 
Fig.Hwe compare and contrast the 7 factors as a function 
of the radial coordinate for the adiabatic and the fully radia- 
tive cases: the dashed (continuous) lines correspond to the 
approximate power-law (exact) solutions. In Fig. [2] we com- 
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FIG. 3. — Comparison between the EQTSs computed using the approximate 
formulas given by Panaitescu & Meszaros 1 1998) (dotted line) and by Sari 
1 1998); Granot et al. 1 1999) (dashed line) in the fully adiabatic case (a = 3/2 
in Eq.|3} and Eq.|S|) and the corresponding ones computed using the exact 
solution of the Taub Eqs.Q (solid line). The difference between the dashed 
line and the dotted line is due to the factor \/2 in the Lorentz 7 factor adopted 
by Sari (see text). The upper (lower) panel corresponds to = 35 s (/„ = 4 
day). 



pare and contrast the EQTS radius at d = as a function of the 
arrival time at the detector. Of particular interest is the cross- 
ing at t~ smaller than 10 2 s between the exact solution and 
the approximate one, which is due to the difference of about 
20 s accumulated over 10 17 cm by the use of the approximate 
Eq.(|6j instead of the exact Eq.(|5}. Note that this difference 
of ~ 20 s in the arrival time, in addition to the effect on the 
EQTSs, impedes the proper identification of the physical pro- 
cesses occurring in the early part of the GRB. 

Two different approximate treatments of the EQTSs exist 
in the current literatur e, bo t h adop ting the above approxima- 
tions a), b), c): one by Sari ( 1998), later used by Granot et al. 

and another one by Panaitescu & Meszaros (1998). 
While lSaril (119981) consid ers only the fully adiabatic case (a = 
3/2 in Eq.Q and Eq.0), Panaitescu & Meszaros ( 1998) con- 
sider both the fully adiabatic (a = 3/2 in Eq.(|3j and Eq.©) 
and th e fully radiative (a = 3 in Eq.(|3j and Eq.0) cases. Sari 
(1998) uses a Lorentz gamma factor T of a shock front prop- 
agating in the expanding pulse, with T = \/2j. 

Without entering into the relative merit of such differing 
approaches, we show in Figs. that both of them lead to 
results different from the ones computed with the exact solu- 
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FIG. 4. — Comparison between the EQTSs computed using the approx- 
imate formulas given by Panaitescu & Meszaros 1 1998 ) (dotted line) in the 
fully radiative case (o = 3 in Eq.l3l and Eq.l6t) and the corresponding ones 
computed using the exact solution of the Taub Eqs.Q (solid line). The upper 
(lower) panel corresponds to tf t = 35 s (t£ = 4 day). 



tions. 

4. COMPARISON WITH THE EXACT SOLUTIONS AND 
CONCLUSIONS 

The consequences of using the approximate formula given 
in Eq.(j3jl to compute the expression t = t(r), instead of the ex- 
act solution of the Taub Eqs.Q, are clearly shown in Figs.|3j- 
|4] The EQTSs represented in these figures are computed at 
selected values of the detector arrival time both in the early 
(~ 35 s) and in the late (~ 4 day) phases of the afterglow. 
Both the fully radiative and fully adiabatic cases are exam- 
ined. Note the approximate expression of the EQTS can only 
be defined for 7 < 7^ and r> r c /. Consequently, at t d a = 35 s the 
approximate EQTSs are represented by arcs, markedly differ- 
ent from the exact solution (see the upper panels of Figs.|3j- 
13. The same conclusion is found for the EQTS at t d = 4 days, 
where marked differences are found both for the fully radia- 
tive and adiabatic regimes (see the lower panels of Figs.|3jS. 

All the observational properties of GRBs, starting from the 
analysis of the prompt radiation ( Ruffi ni et alJl2002l) . to the 
lumin osity in X- and 7- ray bands, to their spectral distribu- 
tion llRuffini et a l. 2003a) as well as inferences on the possible 
presence or absence of beaming in GRBs, depend essentially 
on the structure of the EQTSs. In turn the determination of 
the EQTSs depends on the equations of motion of the bary- 



4 



onic pulse satisfying the Taub equations. The fact that the fi- 
nal results for the observable luminosity, spectral distribution, 
and substructures in the prompt radiation depend on ~ 10 s in- 
tegration paths on different points on the EQTSs implies that 
the agreement between the theoretical predictions and the ob- 
servations becomes a most stringent test for the validity of the 
equations of motion. The correct EQTSs are also essential for 
the identification of th e energy source of the X an d 7 radiation 
in the GRB afterglow (Ruffini etail200l 120031 . 

In conclusion, the approximate treatments largely overesti- 
mate (underestimate) the size of the EQTSs in the early (late) 
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